CARDIVNSWC-TR-94/024  Nonlinear  Calibration  of  an  Infrared  Radiometer 


Carderock  Division,  Navai  Surface  Warfare  Center 

Bethesda,  Maryiand  20084-5000 


CARDIVNSWC-TR-94/024  December  1994 

Signatures  Directorate 
Research  and  Development  Report 


Nonlinear  Calibration  of  an  Infrared  Radiometer 

by 

Peter  O.  Cervenka 

Carderock  Division,  Naval  Surface  Warfare  Center 

and 

Lou  Massa 

City  University  of  New  York 


Presented  at 

Workshop  on  Infrared  Ship  Signature  Modeling  (IRSSM), 
organized  by  the  NATO  Research  Study  Group  on 
Maritime  Infrared  Target  and  Background  Signatures, 
Measurements  and  Characterization  (RSG.5) 

17-18  October  1994 
Paris,  France 


19950321  139 


Approved  for  public  reiease;  distribution  is  uniimited. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OFTHIS  PAGE 


la.REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a.SECURmr  CLASSIFICATION  AUTHORriY 


REPORT  DOCUMENTATION  PAGE 


Ib.RESTRICnVE  MARKINGS 


3.  DISTRIBUnON/AVAILABIUTY  OFREPORT 


2b.  DECLASSIFICATION/DOWNGRADING  SCHEDULE 


Approved  for  public  release;  distribution  is  unlimited. 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER{S) 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


CARDIVNSWC-TR-94/024 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Carderock  Division 
Naval  Surface  Warfare  Center 


6c.  ADDRESS  (Dty,  State,  and  ZIP  Code) 


6b.  OFFICE  SYMBOL 
(If  applicable) 

Code  7230 


7a.  NAME  OF  MONITORING  ORGANIZATION 


7b.  ADDRESS  (CITY,  STATE,  AND  ZIP  CODE) 


Bethesda,  MD  20084-5000 


8a.  NAME  OF  FUNDING/SPONSORING 
CKGANIZATION 


8c.  ADDRESS  (City.  State,  and  ZIP  code) 


6b.  OFFICE  SYMBOL 
(If  applicable) 


11.  TITLE  (Include  Security  Classification) 

Nonlinear  Calibration  of  an  Infrared  Radiometer 


12.  PERSONAL  AUTHOR(S) 

Cervenka,  Peter  O.  and  Massa,  Lou _ ' 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION 

NUMBER 

1  10.  SOURCE  OF  FUNDING  NUMBERS  | 

PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

ACCESSION  NO. 

0602121N 

RH21C17 

3 

13a.  TYPE  OFREPORT 

Final 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED  14.  DATE  OFREPORT  (Year,  Month,  Day)  IS.  PAGE  COUNT 

FROM:  TO:  1994,  December _ 33  _ 


Presented  at  the  Workshop  on  Infrared  Ship  Signature  Modeling  (IRSSM),  organized  by  the 


NATO  Research  Study  Group  on  Maritime  Infrared  Target  and  Background  Signatures,  Measurements  (see  reverse  side) 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identity  by  block  number) 

Dimensionless  Variables,  Infrared  Engineering, 

Infrared  Scaling  Laws,  Scale  Model  Testing,  Thermal  Physics 


19.  ABSTRACT  (Continue  onreverse  if  necessary  and  identify  by  block  number)  | 

At  best,  the  csdibration  of  an  infrared  radiometer,  also  known  as  a  Forward  Looking  Infrared  (FLIR)  system,  is  a  | 
complex  procedure.  Effective  calibration  must  (1)  be  tailored  to  the  actual  radiometer  being  studied,  (2)  be  carried  out  to  each 
of  the  many  subsystems  involved,  (3)  relate  internal  system  parameters  to  physical  quantities  familiar  to  the  infrared 
community,  and  (4)  be  revisited  at  some  interval  to  ensure  that  it  is  still  valid.  The  details  of  some  of  the  ideas  tiiat  form  the 
basis  for  a  valid  calibration  of  the  GEC  Avionics  Ltd  Dual  Waveband  Imaging  Radiometer  FLIR  are  discussed.  The  need  for 
the  removal  of  a  systematic  error  due  to  magnification  optics  emission,  as  described  by  Braim,  is  emphasized.  Another 
systematic  error,  associated  with  an  analysis  of  FUR  gray  scale  response  based  upon  linear  interpolation  firom  reference 
temperatures,  is  examined.  A  novel  method  of  nonlinear  interpolation  is  presented  that  is  aimed  at  the  reduction  or 
elinunation  of  systematic  calibration  errors. 


ZO.DISTRIBUnON/AVAILABllJTY  OFABSTRACT 

UNCLASSIFIEDAJMJMITCD  □  SAMEASRPT. 


22a.  NAMEOF  RESPONSIBLE  INDIVIDUAL 

_ _ Peter  O.  Cervenka  _ 

DDFORM  1473,  84  MAR 


21  .ABSTRAa  SECURITY  CLASSIFICATION 

DTCUSERS  UNCLASSIFIED _ 


22b.  TELEPHONE  (Include  Area  Code) 

f301)  227-1903 


22c.  OFFICE  SYMBOL 

Code  7230 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OFTHIS  PAGE 


UNCLASSIFIED 

SECURmr  CLASSIFICATION  OFTHIS  PAGE 


Block  16  (Continued) 

and  Characterization  (RSG.5),  17-18  October,  1994,  in  Paris,  France. 


I 

DD  FORM  1473,  84  MAR 


UNCLASSIFffiD 


SECURFTY  CLASSIFICATION  OFTHIS  PAGE 


CONTENTS 


Page 

Abstract  .  1 

Administrative  Information .  1 

Introduction .  1 

Calibration  of  GEC  DUWIR  FUR .  1 

FLIR  Description  .  3 

FLIR  Measurement  .  13 

Gray  Scale  Calibration  by  Linear  Interpolation  .  17 

Gray  Scale  Calibration  by  Nonlinear  Interpolation  .  20 

Conclusions  and  Recommendations . .  .  25 

Bibliography .  27 


nCURES 


1.  Typical  calibration  curve  correlating  measured 

values,  y,  to  standard  values,  x .  2 

2.  Measurement  accuracy  and  precision .  4 

3.  FLIR  image  of  a  naval  vessel  at  sea  .  5 

4.  FLIR  system  block  diagram . . .  6 

5.  Single  frame  of  digital  FLIR  data .  7 

6.  FLIR  display .  9 

7.  FLIR  display  and  temperature  data .  10 

8.  Dual  waveband  FLIR  detector .  11 

9.  FLIR  operational  scheme . . .  12 

10.  Internal  reference  surface  temperature  calibration  of  Braim  . . .  18 


11.  Estimated  errors  in  target  apparent  temperature  due  to 
linear  interpolation . 


12.  Comparison  of  various  orders  of  interpolation 


111 


For  j 

mis  GRASil 

DTIC  TAB 

□ 

Unajinotjxiced 

□ 

Justification _ 

1  By - - 

I  Dlgtrlbutlonj^v.  -  _ 

j  Availability  Cbdae 
i  lAvasl  aad/oir 


Special 


8 


21 

24 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


IV 


ABSTRACT 


At  best,  the  calibration  of  an  infrared  radiometer,  also  known  as 
a  Forward  Looking  Infrared  (FLIR)  system,  is  a  complex  procedure. 
Effective  calibration  must  (1)  be  tailored  to  the  actual  radiometer  being 
studied,  (2)  be  carried  out  for  each  of  the  many  subsystems  involved,  (3) 
relate  internal  system  parameters  to  physical  quantities  familiar  to  the 
infrared  community,  and  (4)  be  revisited  at  some  interval  to  ensure  that  it 
is  still  valid.  Vie  details  of  some  of  the  ideas  that  form  the  basis  for  a 
valid  calibration  of  the  GEC  Avionics  Ltd  Dual  Waveband  Imaging 
Radiometer  (DUWIR)  FLIR  are  discussed.  The  need  for  the  removal  of  a 
systematic  error  due  to  magnification  optics  emission,  as  described  by 
Braim,  is  emphasized.  Another  systematic  error,  associated  with  an 
analysis  of  FLIR  gray  scale  response  based  upon  linear  interpolation  from 
reference  temperatures,  is  examined.  A  novel  method  of  nonlinear 
interpolation  is  presented  that  is  aimed  at  the  reduction  or  elimination  of 
systematic  calibration  errors. 


ADMINISTRATIVE  INFORMATION 

This  work  is  in  partial  fulfillment  of  Milestone  2,  Task  3,  of  the  Topside 
Signature  Reduction  Project  (RH21C17)  of  the  Surface  Ship  Technology  Program 
(SC1A/PE0602121N).  The  work  described  herein  was  sponsored  by  the  Office  of  the 
Chief  of  Naval  Research  (OCNR  33)  and  was  performed  by  the  Carderock  Division, 
Naval  Surface  Warfare  Center  (CDNSWC,  Code  7230). 


INTRODUCTION 


The  interest  in  the  FLIR  (Forward  Looking  Infrared)  as  an  imaging  radiometer  is 
focused  in  this  report  on  its  use  in  measuring  the  infrared  contrast  between  a  ship  and  its 
background.  Unless  a  FLIR  is  calibrated,  the  results  obtained  will  be  only  of  a 
qualitative  value.  "The  hundreds  if  not  thousands  of  next  to  useless  uncalibrated 
videotapes  held  across  the  world  containing  valuable  trials  information  whose  data  is  in 
effect  irreducible  bear  witness  to  this.  Calibration  of  the  FLIR  is  therefore  of  crucial 
importance  in  an  IR  ship  measurements  program. 


CALIBRATION  OF  GEC  DUWIR  FLIR 

With  any  calibration,  the  objective  is  to  obtain  a  correlation  of  measured  values, 
y,  with  true  (standard)  values,  x.  Typically,  such  a  correlation  is  represented  by  a 
"calibration  curve"  such  as  that  shown  in  Fig.  1.  The  possibility  of  calibration  is  based 
upon  comparison  to  a  trusted  standard.  Confidence  in  a  calibration  is  limited  by  that  of 
the  accessible  standard. 
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Fig.  1.  Typical  calibration  curve  correlating  measured  values,  y, 
to  standard  values,  x. 
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Measurements  in  general  and  calibrations  in  particular  may  contain  systematic  and 
random  errors.  As  far  as  possible,  calibration  strives  to  eliminate  systematic  errors  and 
to  place  an  upper  boundary  on  random  errors.  Systematic  errors  (sometimes  due  to 
instrument  deficiency)  and  random  errors  (for  instance  due  to  variable  fluctuations)  affect 
measurement  accuracy  and  precision,  respectively.  Accuracy  is  measured  by  the 
"distance"  between  a  test  estimate  of  a  physical  variable,  such  as  the  average,  and  its 
standard  value;  see  Fig.  2. 

One  purpose  of  calibration  is  to  provide  an  understanding  of  the  physical  sources 
of  measurement  errors.  For  this,  it  is  required  to  create  a  theory  that  encompasses  both 
the  variable  measured  as  well  as  the  instrumentation  being  used.  If  at  all  possible,  the 
aim  is  to  reduce  both  error  types  and  to  estimate  the  magnitudes  of  the  errors  that  remain. 
Such  understanding  can  then  form  the  basis  for  an  effective  use  of  the  measured  data  and 
for  subsequent  improvement  in  the  calibration. 

Unfortunately,  calibrations  are  never  concluded  once  and  for  all.  Instrument 
characteristics  are  altered  by  time  and  use,  and  calibration  must  be  viewed  as  an  iterative 
process  as  long  as  an  instrument  is  in  use.  It  is  good  practice  for  a  measurement 
laboratory  to  maintain  reliable  standards  against  which  periodic  recalibration  may  be 
performed. 


FUR  DESCRIPTION 

The  Carderock  Division  of  the  Naval  Surface  Warfare  Center  in  Bethesda, 
Maryland,  operates  a  Dual  Waveband  Imaging  Radiometer  (DUWIR)  manufactured  in  the 
United  Kingdom  by  GEC  Avionics  Ltd.  The  following  analysis  is  applicable  to  the 
DUWIR  system. 

The  measurement  produced  by  a  FLIR  is  a  thermal  map  of  an  infrared  radiation 
source  and  its  surroundings.  Calibration  will  affect  the  entire  image  produced  by  a 
radiometer.  Figure  3  shows  such  a  thermal  scene  in  the  8-  to  12-/im  infrared  band 
(LWIR),  which  represents  a  naval  vessel  at  sea  on  a  sunny  day. 

During  the  measurement  process,  imagery  is  first  recorded  as  analog  video  on 
magnetic  tape.  Figure  4  shows  that  the  FLIR  system  also  contains  image  analysis 
hardware  that  allows  the  capture  in  digital  format  of  any  single  image  frame.  A  typical 
single  frame  of  FLIR  data  in  digital  format  is  shown  in  Fig.  5.  The  image  is  built  up 
from  a  rectangular  array  of  pixels.  The  signal  voltage  from  the  FLIR  for  each  pixel  is 
represented  by  an  8-bit  number  capable  of  being  represented  on  a  gray  scale  G  of  256 
discrete  levels.  The  intensity  or  gray  scale  level  of  each  pixel  may  be  examined 
individually,  and  associated  with  a  temperature  and/or  radiance  value  for  that  particular 
pixel.  Other  data  include  two  reference  internal  temperatures,  displayed  in  Fig.  5  as  Pelt. 
1  and  Pelt.  2,  as  well  as  a  scanner  internal  temperature.  In  addition,  the  width  of  the 
temperature  window  W  and  the  offset  O  which  sets  the  temperature  level  are  also  output. 
These  and  other  data  are  displayed  around  the  FLIR  image  in  Figs.  6  and  7. 
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Fig.  4.  FLIR  system  block  diagram 
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Fig.  5.  Single  frame  of  digital  FLIR  data. 
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The  GEC  Avionics  Ltd  FLIR  uses  an  infrared  detector  (Fig.  8)  that  includes  two 
SPRITE  (Signal  Processing  Right  in  the  Element)  components  respectively  sensitive  in  the 
4-  to  5-/rm  and  8-  to  12-/xm  wavebands.  The  detector  material  for  both  wavebands  is 
mercury  cadmium  telluride. 

The  FLIR  operational  scheme  is  shown  in  Fig.  9.  A  study  of  this  block  diagram  is 
especially  useful  in  understanding  FLIR  calibration.  Note  that  the  image  radiation  and 
optical  path  are  identical  for  both  wavebands  from  target  to  detector.  Beyond  the 
detector,  separate  electronic  amplification  and  signal  processing  serve  to  record  and 
display  the  images  in  each  waveband.  Note  also  that  two  Peltier  blackbody  reference 
surfaces  are  viewed  at  the  edges  of  each  image  frame.  Platinum,  Pt,  resistance  devices 
measure  the  temperatures,  Tj  and  T2,  of  the  reference  surfaces.  Not  shown  in  the  figure 
is  a  Pt  resistance  derivative  which  monitors  the  temperature,  Tg,  of  the  FLIR  scanning 
components.  All  of  these  temperatures  displayed  in  Figs.  3  and  7  through  9  are  involved 
in  the  FLIR  calibration.  In  fact,  FLIR  calibration  minimally  involves  the  measurement 
variables  T^,  T2,  and  T^;  that  is,  the  gray  scale  response  of  the  pixel  (m  and  n)  for 
all  m  and  n,  the  temperatures  of  reference  surfaces  1  and  2,  and  the  temperature  of  the 
FLIR  scanner.  As  indicated  in  Fig.  5,  each  pixel  gray  scale  response  Gj^j,  is  to  be 
associated  with  a  temperature,  T^„.  It  is  possible  to  consider  that  the  basic  variable  of  the 
FLIR  calibration  is  the  temperature  or,  more  specifically,  several  temperatures,  Tj„„  for 
all  (m  and  n),  T^,  T2,  and  T^.  The  goal  in  these  measurements  is  to  eliminate  systematic 
errors  and  to  quantify  random  errors. 

Obtaining  three  of  the  quantities, Tj,  T2,  and  T^,  is  fairly  straightforward  since  the 
procedure  in  each  case  involves  a  contact  measurement  of  temperature.  The  calibration 
simply  requires  correlation  with  a  known  standard,  such  as  a  standard  thermocouple  put 
into  contact  and  thermal  equilibrium  with  the  same  thermal  mass  as  the  one  being 
measured.  In  the  case  of  Tj  and  T2,  the  thermal  mass  is  that  of  the  internal  reference 
surfaces  in  the  FLIR;  and  in  the  case  of  T^,  that  of  the  FLIR  scanner  assembly.  In  each 
of  these  cases  the  calibration  is  a  simple  procedure,  resulting  in  a  single  calibration  curve 
for  each  of  the  three  temperature  measurements,  similar  to  that  shown  in  Fig.  1. 

The  temperature  of  the  two  internal  reference  surfaces,  Tj  and  T2,  should  be 
calibrated  in  two  different  ways.  One  measurement  is  a  contact  measurement  which  gives 
the  equilibrium  temperature  of  the  bulk  mass  of  the  internal  reference.  It  is  this 
measurement  that  would  be  calibrated  against  a  standard  thermocouple  attached  to  the 
same  thermal  mass.  Another  measurement  incorporates  the  gray  scale  response  of  the 
FLIR  viewing  an  external  radiating  standard  blackbody  of  known  temperature.  The  gray 
scale  response  of  the  FLIR  to  the  external  standard  is  matched  against  that  of  the  internal 
reference  to  obtain  a  surface  temperature  of  the  internal  reference.  This  surface 
temperature  calibration  is  more  complex  than  the  surface  contact  measurement. 
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Fig.  8.  Dual  waveband  FLIR  detector. 
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Fig.  9.  FLIR  operational  scheme.^ 
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FLIR  MEASUREMENT 


To  effectively  calibrate  the  FLIR,  a  rudimentary  theory  must  be  constructed  that 
includes  the  FLIR  measurement  as  well  as  its  errors.  Measurements  must  then  be  carried 
out  to  test  the  theory  and  to  assign  magnitudes  to  the  errors.  If  the  errors  given  by  theory 
and  experiment  are  about  the  same  magnitude  and  follow  similar  trends  as  a  function  of 
experimental  variables,  the  errors  may  be  considered  as  being  reasonably  well  understood 
and  accounted  for. 

Fortunately,  a  very  good  theory  in  first  approximation  has  been  constructed  for 
the  FLIR  by  Braim.^  In  his  treatment,  the  FLIR  is  assumed  to  view  a  large  target  which 
radiates  in  accordance  with  Planck’s  blackbody  formula: 

C\  /i\ 

W(T)  =  — - 1 -  (1) 

X^[exp(C2/Xr)-l] 


where: 


W 

T 

X 

C 


1 


emitted  flux,  W/m^//xm 
absolute  temperature,  K 
wavelength,  fim 
3.7415  X  10®,  urn^ 
1.4388  X  10^  fimK 


The  FLIR  detector  produces  a  signal  voltage,  S,  which  depends  upon  the  received 
flux,  <i),  according  to 


5= 


(2) 


where  R  is  the  detector  responsivity.  The  flux  that  reaches  the  detector  from  the  target 
has  contributions  from  all  components  in  the  optical  path  originating  at  the  target,  t. 
These  components  are:  the  atmosphere,  a;  the  magnification  optics,  o;  and  the  FLIR 
scanner,  s.  Each  of  these  components  is  assumed  to  have  an  emissivity 


(3) 


where  a  is  the  component  absorptivity.  If  for  each  component. 


a+T=l  , 


(4) 


where  t  is  the  component  transmissivity,  then 


6=1 -r  .  (5) 

If  this  expression  holds  for  the  emissivity  of  each  component,  then  the  flux  generated  by 
each  component  along  the  optical  path  is 


(6) 

4>0=<t>a^O+MTo)(l-Ta)  , 

(7) 

<^,=V,+W(7;)(1-7,)  . 

(8) 

Insertion  of  these  component  contributions  to  the  flux  into  Eq.  (2)  results  in  a  signal 
voltage  at  the  detector  due  to  the  target 


5,=  j  /?W(r,)T„ToT^+ 1  W(7;)(l  -tJtqt/K 

+  1  i?W(ro)(l  /?W(r,)(l  -T^)dK 


In  this  equation  the  four  integrals  on  the  right  hand  side  represent  emission  from  the 
target,  the  atmosphere,  the  magnification  optics,  and  the  scanner,  respectively.  The  flux 
that  reaches  the  detector  from  the  internal  reference  surface  contains  contributions  from 
the  components  in  the  optical  path  originating  at  the  reference  surface.  In  this  case  the 
only  such  component  is  the  FLIR  scanner.  Beyond  this  component,  the  flux  is 


(10) 


Inserting  this  expression  into  Eq.  (2)  yields  the  signal  voltage  at  the  detector  due  to  the 
reference  surface 


S,^lRWa,^r/!k*^RW{T,)(l-7,)dK  .  (11) 

The  two  integrals  on  the  right-hand  side  of  Eq.  (11)  represent  the  emission  from  the 
reference  surface  and  from  the  scanner,  respectively.  A  temperature  reading  given  by  a 
platinum  wire  resistance  in  contact  with  a  reference  surface  is  obtained  by  adjusting  the 
Peltier  device  of  the  surface  so  that  the  gray  scale  response  of  the  target  and  reference 
surface  match.  Thus, 
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(12) 


Rearranging  slightly  and  using  Eqs.  (9)  and  (11) 

I  RW(r,)T^ToT/x= I  I  i?w(r„)(i  -t^)to7^ 

-|RW(ro)(l-To)rA 

Note  that  under  conditions  of  the  same  gray  scale  response  from  target  and  internal 
reference  surface,  the  emission  from  the  target  is  equal  to  that  from  the  reference  surface 
reduced  by  contributions  from  the  atmosphere  and  the  magnification  optics.  In  case  Eq. 
(12)  holds,  there  is  no  contribution  due  to  scanner  emission. 

Equation  (13)  may  be  used  for  FUR  calibration.  Imagine  for  simplicity  an  ideal 
case  whereby 


(14) 


In  this  case,  Eq.  (13)  reduces  to 

I  RW(T,)t/!K^  I  RWiT,^)T^  .  (15) 

Thus,  if  the  atmosphere  and  optical  train  do  not  result  in  any  detectable  emissions, 
matching  the  gray  scale  responses  of  target  and  reference  surface  also  matches  their 
respective  temperatures.  This  is  a  case  where  it  is  conceptually  simple  to  calibrate  T^^f. 
For  a  range  of  target  temperatures  Tj,  the  reference  surface  Peltier  device  is  adjusted  so 
that  target  and  reference  surfaces  both  have  the  same  gray  scale  responses.  Then 
Tt  =  T^f 

On  the  other  hand,  if 


(16) 


then  Eq.  (13)  reduces  to 


I  RW(r,)ToT^=|  RW(T,^)t/K 
-|RW(ro)(l-To)T^ 


(17) 
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The  last  integral  on  the  right-hand  side  of  Eq.  (17)  represents  a  correction  term  to 
the  first  integral  on  the  same  side;  that  is,  it  plays  the  role  of  a  systematic  error. 

Matching  the  gray  scale  response  of  target  and  reference  surface  no  longer  implies  that 
Tj  =  Tjgf ,  as  was  the  case  previously.  Although  this  would  still  be  approximately  true, 
the  actual  temperature  ,  for  given  values  of  Tj  and  T^,  can  be  obtained  by  using  an 
iterative  solution  to  Eq.  (17). 

Different  values  of  T^,  will  modify  the  solution  obtained  for  T^gf,  signifying  that 
the  temperature  of  the  magnification  optics  introduces  a  systematic  error  into  the 
determination  of  by  the  process  of  requiring  Eq.  (12)  to  hold.  The  magnitude  of  that 
systematic  error  is  controlled  by  the  last  integral  in  Eq.  (17).  As  the  transmissivity  of  the 
optical  train  approaches  unity  (that  is,  as  the  emissivity  vanishes),  so  too  does  the 
systematic  error,  and  Eq.  (17)  reverts  back  to  Eq.  (15). 

Considering  Eq.  (13)  in  light  of  these  last  remarks,  it  appears  that  the  integrals 
subtracted  on  the  right-hand  side  describe  the  atmospheric  emission  and  the  optical 
emission,  respectively.  These  are  the  systematic  errors  due  to  temperature  matching  of  the 
gray  scale  response  of  the  target  and  the  internal  reference  surface,  respectively.  To  show 
in  a  quantitative  manner  that  Eq.  (13)  represents  a  valid  description  of  the  temperature 
measurement  and  its  errors,  requires  that  the  quantities  in  question  be  both  calculated  and 
measured.  Braim  has  done  that  for  the  8-  to  12-;im  waveband;  see  Fig.  10.  Figure  10a 
shows  the  experimentally  measured  results  of  T^gf  versus  Tj,  for  a  variety  of  different 
temperatures,  T^  of  the  optical  train.  Note  that  with  increasing  T^,  T^gf  is  forced  to 
increase  to  compensate  for  the  effect  of  increased  radiation  from  the  optical  train,  in 
accordance  with  Eq.  (13).  Part  (B)  of  the  figure  presents  theoretically  calculated  results, 
based  on  Eq.  (13),  of  T„f  versus  T^,  for  a  variety  of  different  optical  train  temperatures, 
Tq.  a  clear  correlation  exists  between  the  experimental  and  theoretical  results,  which 
lends  confidence  to  Eq.  (13)  as  a  good  first-order  description  of  the  FLIR  temperature 
measurement  and  its  errors. 

Although  the  full  calibration  described  by  Eq.  (13)  takes  into  account  the 
atmospheric  transmissivity,  t^,  and  temperature,  Tg,  these  quantities  are  not  considered  in 
Fig.  10.  This  is  a  matter  of  first-order  expediency,  and  a  more  complete  treatment  should 
also  take  into  account  systematic  errors  due  to  the  atmosphere. 

Greater  understanding  of  the  FLIR  calibration  process  has  been  gained  through 
measurement  of  T^f.  The  FLIR  produces  a  gray  scale  response  to  a  blackbody  target  with 
known  temperature  which  is  considered  to  be  a  standard  value.  The  temperature  of  the 
Peltier  internal  reference  surface  is  adjusted  so  that  its  gray  scale  response  matches  the  , 
standard,  a  condition  that  is  embodied  in  Eq.  (13).  The  resistance  measurement  of  the 
temperature,  T^f,  of  the  internal  reference  via  a  platinum  wire  resistance,  is  correlated  to 
the  blackbody  target  temperature,  \  This  correlation  is  established  over  a  temperature 
range  of  interest  for  a  fixed  value  of  the  magnification  optical  components  temperature, 

T  .  In  order  to  remove  the  systematic  error  associated  with  the  optical  train  temperature, 
T°,  in  Eq.  (13),  T^gf  versus  Tj  is  measured  for  a  range  of  temperatures,  T^,  building  a  set 
of  experimental  calibration  curves  similar  to  those  displayed  in  Fig.  10a.  Such  curves 
allow  the  measurement  of  an  unknown  apparent  temperature,  Tj,  of  a  target  by  means  of 
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matching  the  internal  reference  temperature  gray  scale  responses.  It  may  thus  be  possible 
to  remove  the  systematic  error  due  to  T^. 

Calibration  of  the  internal  reference  surface  temperature,  Tj.gf,  is  adequate  for 
measuring  the  temperature,  Tj,  of  an  extended  target  with  uniform  temperature.  Many 
targets  of  interest  are  characterized  by  a  nonuniform  temperature  distribution,  which 
corresponds  to  a  number  of  different  gray  scale  responses.  This  case  requires  the 
calibration  of  the  entire  gray  scale  for  every  pixel  of  an  image  frame.  The  following  gray 
scale  calibration  is  based  upon  the  calibration  of  both  internal  reference  surfaces  using  the 
method  already  presented. 


GRAY  SCALE  CALIBRATION  BY  LINEAR  INTERPOLATION 

Assume  that  the  temperatures  of  the  two  internal  reference  surfaces  are  calibrated 
using  the  theory  shown  in  Fig.  10  and  consistent  with  Eq.  (13).  For  each  discrete 
temperature  of  a  reference  surface  there  is  a  corresponding  gray  scale  response  curve.  For 
each  set  of  reference  surface  temperatures  there  is  a  corresponding  gray  scale  response 
curve.  The  desire  is  therefore  for  an  interpolation  scheme  that  assigns  temperatures  to  the 
gray  scale  data  of  a  given  image  frame. 

Assuming  that  the  gray  scale  response,  G,  is  a  function  of  temperature,  T,  which 
to  first-order  may  be  expanded  by  a  Taylor-series  of  the  form 


G{T)  =  (7(ri)+(7’-ri)_^  ,  (18) 

where  T^  and  T2  are  the  temperatures  of  the  two  reference  surfaces.  A  finite  difference 
evaluation  of  the  first  derivative  of  Eq.  (18)  produces 


bG{T{)  _  G{T2)-G{T{) 
5T  T2-Ti 


The  Taylor-series  expansion  becomes 


G{T)=G(T^)+{T-T^) 


G{T2)-G{T{) 

T2-T1 


(20) 


which,  if  solved  for  T,  yields 
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Referonce  tefnp«rotur«  (*C) 


Fig.  10a.  Experimental  measurements. 


Fig.  10b.  Corresponding  theoretical  results. 


Fig.  10.  Internal  reference  surface  temperature  calibration  of  Braim. 
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(21) 


Defining 


T=T^^[G{T)-G{T^)] 


(7’2-7’i) 

[G{T2)-g{T{)]  ■ 


M=[G{T)-G{T^)] 


ih-T^) 

{G{T2)-G{T{)] 


(22) 


obtains 


T=Ty+LT  .  (23) 

The  gray  scale  response,  G,  of  a  given  pixel  in  an  image  may  thus  be  used  to  obtain  the 
temperature,  T,  of  that  same  pixel. 

Since  the  Planck  blackbody  curve  is  not  linear  in  temperature,  a  linear 
interpolation  scheme  will  increasingly  become  more  unreliable  as  the  separation  of  the 
two  reference  temperatures,  T2  and  T^,  increases.  The  magnitude  of  the  error  due  to 
linear  interpolation  may  be  estimated  by  calculating  the  signal  voltage 


S,=|i?W(rpi/X  /=1,2  (24) 

at  the  two  reference  temperatures,  and  then  taking  the  "midpoint"  signal 

S=1(S,+S2)  ,  (25) 

which  occurs  at  the  midpoint  temperature 

r=l(ri+72)  •  (26) 


S  may  be  obtained  through  the  use  of  a  relationship  between  the  signal  voltage  and  the 

temperature,  such  as  Eq.  (24).  The  results  of  the  calculations  performed  in  this  manner  by 
Braim  are  shown  in  Fig.  11.  The  error  due  to  linear  interpolation  of  the  gray  scale  is 
seen  to  increase  with  the  temperature  interval,  T2  -  Tj,  and  can  easily  amount  to  a  few 
degrees  absolute.  In  the  shorter  waveband,  the  corresponding  errors  are  noticeably  larger. 
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GRAY  SCALE  CALIBRATION  BY  NONLINEAR  INTERPOLATION 

The  use  of  a  linear  interpolation  scheme  introduces  an  extraneous  systematic  error 
into  the  final  result.  More  data  exist  in  the  gray  scale  than  are  used  by  a  linear 
interpolation  method.  A  linear  interpolation  scheme  establishes  a  linear  relationship 
between  the  signal  received  by  the  detector  from  an  element  of  the  target  and  the 
temperature  of  that  target  element.  The  detector  responds  to  radiation  described  by  the 
Planck  blackbody  expression  which  is  not  linear  in  T.  Thus  a  systematic  error  is 
introduced  by  the  linear  interpolation  scheme.  This  systematic  error  can  be  eliminated  or 
reduced  through  the  use  of  nonlinear  interpolation. 

The  platinum  resistance  temperature  reading  of  an  internal  reference  surface,  of 
course,  is  independent  of  the  changing  electronic  "window"  and  "offset'  values.  The 
reference  surface  temperatures  are  thus  calibrated  against  a  blackbody  reference  source. 
Given  such  a  calibration,  the  task  now  is  to  use  that  information  to  calibrate  the  gray 
scale  without  introducing  the  systematic  error  associated  with  a  linear  interpolation. 

The  gray  scale  response  of  a  reference  surface  ideally  is  linear  in  the  signal 
voltage  induced  in  the  detector  by  radiation  from  a  surface;  that  is. 


where  k  is  a  proportionality  constant,  the  magnitude  of  which  is  fixed  by  the  amplification 
electronics  of  the  FLIR,  W  is  the  Planck  blackbody  function,  R  is  the  response  function 
of  the  detector,  Tj  is  a  reference  surface  temperature,  and  Tg  is  the  temperature  of  the 
scanner.  The  constant  satisfies 


k  =  kiW,0)  ,  (28) 

that  is,  the  constant  is  a  function  of  the  electronic  "window"  W  and  "offset"  O 
parameters.  This  makes  the  gray  scale  response  inconvenient  to  calibrate  directly  against  a 
blackbody  source.  For  any  given  frame  of  image  data  the  constant  k  may  be  evaluated 
from  Eq.  (27).  Thus, 


kiW,0)= - - -  .  (29) 

|R[W(r,)+W(r5)(l-T5)]t»v 

Such  determination  of  k  requires  that  the  detector  responsivity  and  the  scanner 
transmissivity  have  been  measured.  From  the  Eq.  (29),  the  value  of  k  can  be  used  to 
solve  for  the  temperature,  T,  to  be  assigned  to  every  pixel  in  the  image  frame.  For  any 
pixel  the  gray  scale  response  is 
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TEMPERATURE  (K)  TEMPERATURE  <K) 

Fig.  11a.  5-nm  wide  spectral  window.  Fig.  11b.  12-/im  wide  spectral  window. 


Fig.  11. 


Estimated  errors  in  target  apparent  temperature 
due  to  linear  interpolation.^ 
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(30) 


G(7)=it|/?[W(7)+W(r5)(l-T5)]^/X  . 


This  equation  may  be  solved  iteratively  for  the  temperature  of  the  pixel,  circumventing 
any  linear  interpolation  scheme. 

An  iterative  solution  to  this  equation  may  be  considered  somewhat  inconvenient. 
In  which  case,  an  interpolation  scheme,  which  include  the  information  contained  in  the 
constant  k  and  given  by  Eq.  (29),  would  be  somewhat  superior  to  a  linear  one. 

Instead  of  using  Eq.  (30)  to  solve  for  all  the  values  of  T  and  all  pixel  values  G,  it 
is  used  only  once  to  solve  for  G  corresponding  to  a  T  value  half  way  between  two 
reference  surface  temperatures,  Tj  and  T2.  Thus, 


(31) 


where 


h  =  icJ-j-r,) 


(32) 


If  one  Taylor-expands  the  gray  scale  response,  G,  to  second-order. 


G(7)=r(7’i)+(r-ri) 


dGiT^)  ^  (T-T^f  d^GiT^) 
6T  2  572 


(33) 


Inserting  T,,  into  Eq.  (30)  results  in  G  for  three  temperatures,  Tj,  T^,,  and  T2.  Knowledge 
of  G  at  three  temperatures  allows  an  evaluation  of  both  the  first  and  second  derivatives  by 
finite  difference  approximations  given  by 


G{T,^h)-GiT,) 

h 


(34) 


and 


d^GiT^) 

6T^ 


GiTi+2h)  -  2GiTi+h)+G(T{) 


(35) 


respectively. 
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If  the  gray  scale  response,  G,  is  Taylor-series  expanded  to  second-order,  then 


GO)  =  C7(ri)+(r-ri) 


bG{T{)  ^  b^G{T{) 

5T  2  572 


(36) 


using 


AT=(T-Ti) 


(37) 


The  Taylor-series  expansion  is  a  quadratic  equation  in  the  variable  AT;  that  is. 


(A7)2 


1  8^G 

2  8T 


(^i) 


+Ar 


6G(ri) 

8T 


-[GiT^)-G(T)]=0 


(38) 


with  solution 


AT=. 


-8G 

8T 


(7’i)± 


— (^i) 

8T  ^ 


16W1) 
2  572 


[G(ri)-G(7)] 


Vi 


6^G(ri) 

8T^ 


(39) 


Using  Eq.  (39),  the  temperature  of  any  pixel  with  gray  scale  response  G(T),  becomes 


T=Ti+AT  . 


(40) 


Equation  (40)  is  a  nonlinear  interpolation  formula  valid  to  second-order  in  the  Taylor- 
series  expansion  of  the  gray  scale  response. 

Figure  12  indicates  qualitatively  the  nature  of  the  improvement  expected  when 
using  the  nonlinear  second-order  interpolation  formula,  Eqs.  (39)  and  (40X  instead  of  the 
linear  first-order  interpolation  formula,  Eqs.  (22)  and  (23).  The  points  in  Fig.  (12), 
labeled  1,  2,  and  3  refer  to  a  Taylor-series  expansion  evaluated  up  to  order  zero,  one, 
and  two,  respectively.  Examination  also  indicates  two  reference  data  points  determine  the 
first  derivative  of  a  function  and  three  reference  data  points  determine  the  second 
derivative,  respectively. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

A  FLIR  can  be  used  as  an  imaging  radiometer  for  quantifying  the  infrared  contrast 
between  a  ship  and  its  background  will  have  qualitative  value  only  if  care  is  taken  in  the 
calibration  of  the  instrument. 

Following  the  theory  of  Braim,  the  importance  of  the  removal  of  one  systematic 
measurement  error  due  to  emission  from  the  FLIR  optical  train  is  explained.  Another 
source  of  systematic  error,  attributable  to  analysis  of  the  FLIR  gray  scale  response  and 
based  upon  a  method  of  linear  interpolation,  may  be  eliminated. 

A  novel  analysis,  both  nonlinear  and  iterative,  using  a  second-order  Taylor  series 
expansion,  is  presented.  It  is  believed  that  this  method  may,  at  least,  largely  reduce  the 
calibration  systematic  error  introduced  by  the  linear  interpolation  method  of  calibration. 

Further  analysis  is  required  to  analyze  the  nature  of  the  random  errors  associated 
with  FLIR  measurements. 
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